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Abstract 

We study the radiative pion decay of vr^ — )• e^VeJ in the hght front quark model (LFQM). We 
also summarize the result in the chiral perturbation theory. The vector and axial-vector hadronic 
form factors (Fy^j^) for the vr — )• 7 transition are evaluated in the whole allowed momentum transfer. 
In terms of these momentum dependent form factors, we calculate the decay branching ratio and 
compare our results with the experimental data and other theoretical predictions in the literature. 
We also constrain the possible size of the tensor interaction in the LFQM. 
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I. INTRODUCTION 



The light pesudoscalar decays have been playing important roles of understanding the 
standard model (SM). In particular, the radiative pion decay of 7r+ — )■ e^Ugj {'iTe2-y) is an 
interesting process, which can be used to test the V ~ A structure of the weak interaction 
and search for some anomalous interactions beyond the SM. The decay consists of two 
types of contributions, referred as internal-bremsstrahlung (IB) and structure-dependent 
(SD) in terms of the emission of the photons, respectively. The IB contribution to the 
decay amplitude [Mib) is helicity suppressed like the 'Ke2 decay as the photon radiates from 
the external electron, while the SD one [Msd)-, depending on vector and axial- vector weak 
hadronic currents, is proportional to the electromagnetic coupling constant a but free of the 
helicity suppression. One can parametrize Msd by the vector and axial- vector form factors, 
denote as Fy and Fa-, respectively. 

The decay of vr"'' — )■ e^VfTi has been measured with the branching ratio of (1.61± 
0.23) X 10^^ for the cuts of > 21 MeV and E^ > 70 - 0.8E^ by the ISTRA experi- 



ment 



2|]. Recently, a more precise measurement on the decay branching ratio has been 
given by the PIBETA Collaboration js, 4|, with the decay branching ratios in various kine- 
matic regions. In particular, for the cuts of E^. > 0.5 MeV and E^ > 10 MeV with the 
relative angle Oe^y > 40°, the decay branching ratio is (73.86 ± 0.54) x 10^*^ [4]. The new 
ongoing PEN experiment at PSI will at least double the PIBETA data set A -ulting in 
further improvements in precision [6[. In addition, there is another ongoing new experiment, 
PIENU, at TRIUMF j3] with a similar sensitivity as the PEN experiment. 

Theoretical calculations on Fv,a as well as the decay branching ratio in the SM have 



been done in various QCD models 



ll|-l3]. In particular, the decay branching ratio with 



the same cuts as those by ISTRA and PIBETA ^ is found to be 2.55x10 and 

76.66x10^^ in the chiral perturbation theory (ChPT) at 0{p^) 8l-ll0|. which are larger the 
data shown above, respectively. As a result, it may be necessary to consider some new types 
of interactions, such as tensor interactions It is clear that these tensor interactions 

are undoubtedly signals of new physics. On the other hand, it is important if we can obtain 
information on Fy^A in some QCD models other than the ChPT. For this purpose, in this 
study we will evaluate Fy^A in the light front quark model (LFQM) 18|, ll9| . We will use the 



form factors in both ChPT and LFQM to examine the decay of vr^ — t- e"^z/e7. In addition, 



we will examine the new physics effect due to the tensor interactions. 

This paper is organized as follows. In Sec. II, we summarize the form factors in the vr — )■ 7 
transition within the ChPT and LFQM. In Sec. Ill, we calculate the decay branching ratio 
of vr"*" — )• e"'"z/e7 in these models. We also compare our results with the experimental data 
and other theoretical predictions in the literature. We give our conclusions in Sec. IV. 



II. THE FORM FACTORS 



A. Vector and Axial-vector Form Factors 



The decay amplitude for — )■ e"*"z/e7 can be written as: [20|, [21 1 



M = MiB + MsD 



MiB = ie^Vudf-Krriee* u{pe){l - 75) ( — 
V2 KPtt ■ 



Gf 

MsD = -2--^Kde*M(pe)7a(l " l5)v{p^ 



P^ 2p^+ )^7^ 

q 2pe-q 
Fa 



v{Pv 



e— ■ q + Kg" + ^e—e^'-^-qpp^x 



,(1) 



where is the photon polarization vector, Pn, Pe, Pu, and q are the four momenta of 7r+, e"*", 
z/ and 7, and /^r and FAy are the tt meson decay constant and the axial-vector and vector 
form factors, defined by 



{0\sYl5u\7r+{Pn)) =ifnP':,, 

p 

(7(?)|w7/.75f^|7r(p^) ) = e— [{p ■ q)e*^ - {e* ■ p)q^] 



(2) 



respectively, with p = Pn — q being the transfer momentum. Obviously, Mjb has a suppres- 
sion factor of mg. The physically accessible kinematics region is < < p^ax — ''^^ due to 
the time-like momentum transfers. In the following discussion, we will first summarize the 
formulas for FyA in the ChPT BQandthe^ evaluate t.e.efo.. facto. „U.e LFQM. we 



note that similar calculations 
have been performed in Refs. 



br the P ^ 7 (P = K+, D, B) transitions in the LFQM 
22|-3. 



1. Chiral Perturbation Theory 



The tree and loop contributions to Fv,a in the ChPT at 0{p^) for the tTpo ^t- decay have 
been calculated in Refs. jo], [lo|. The explicit forms can be summarized as 22| 
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m 



2V2 



6F3(27r) 
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+ 2m^(62/[7 + Syls - 2|/8i - 2^82 + 22/83 + 1/84 + 2/85 - yloo + 22/[o2 
+ 22/[o3 - 22/U + l/Iog) + P^K^Vwo ' %[o9 + vlio] 



(4) 



where the wave function and decay constant (F^r = /tt/v^) renormalizations have been in- 
cluded and C^, LI and yl are the renormalized coupliiig constants. Note that the first terms 



in Eqs. ([3]) and ([4]) correspond to Fy and Fa at O(p^) [8, 25], respectively. To get the numer 



ical results for the form factors, we take rriK = 0.495 GeV, = 0.14 GeV and rup = 0.77 
GeV, = 0.092 GeV and the renormalized coefficients of (L ^, Lp , L3, Lg, L\q), C22) and 
IZlpo, l/[o4, ylm^ 2/no) ^ be (0.53, 0.71, -2.72, 6.9,-5.5) x lO'^ y, (0.013, 6.52) x 10"^ GeV-^ 



271 and (1.09, -0.36, 0.40, -0.52) x 10-^/^2 |28|, respectively. For some other possible sets 



of coefficients, see Ref. [101 as well as the recent review in Ref. 



29|. Note that the uncer- 



tainties for the renormalized coupling constants are not considered in this study. 



2. Light Front Quark Model 



' the phenomenological LF meson 
22|. The pion wave function can 



In the hght front (LF) approach, the general structure o 
wave function is based only on the Qq Fock space sector 
be expressed by an anti-quark q and a quark Q with the total momentum (p + q) as: 

Hp + (1)) = Y1 J [^^i] [dh]2{2TTf6^{p + q-h-k2) 

XIX2 

X ^t^'^{z,k^)b+{h,X,)d^{h,X2)\0) , (5) 

where jg ^j^g amplitude of the corresponding q{Q) and ki(^2) is the on-mass shell LF 

momentum of the internal quark. The LF relative momentum variables {z, k±) are defined 
by 

kf = {I- z){p + qy, kt = z{p + qy, 

ki± = (1- z){p + q)± + k±, k2± = z{p + q)±- k±, (6) 



and 



/ k+k+ \ ' 

<i>t^'^{z,k^)= / ' ^ u{k^,\l)l'v{k2,\2)H^,k^), 

\2[M^ - {niQ - mg) ]J 

= ^ + — ^ . 7 

1 — z z 

where (f){z, k±) is the space part of the wave function, which is taken to be a Gaussian type 
but it can be solved in principle by the LF QCD bound state equation [24]. At the quark 
loop diagram, the hadronic matrix elements in Eq. (|2]) can be obtained to be 

/d^k 

r i{-^2 + mu). i{I^3 + mu) ,^ i^kfi + md) 
^1^5 ,0 o , .2 o , , 7m(1 - 75) ,2 o 



' kl — m\ + ie kl — m\ + ie kf — 7712 + ie 

+ {u ^ d,ki^ k2)> , (8) 



where is a vertex function related to the quark- ant iquark bound state of the vr meson, 
k2 = q — ks and ki = {p + q) — k2 = k^ + p. By integrating over the LF momentum A;^ in 
Eq. §1, we get 

{l{q)\ul^^ (1 - 75) d\ 7r(p + q) ) 

'^\d'kJ-^nj—i-^ + iu^d,k,^k2)] , (9) 



where 



{ri^u 1 _ dkfdki± 

^ ~ 2{27rrktk+kt ' 

^''Ikz = Trl 75(-^2 + mu)iegf{l^3 + m„)7^(l - 75)(^i + ma) 



kion — ; ^1(2) ~ Pan ~ ^2(l)on ' ^3 ~ 9' ^ ^lon i (lO) 

with {on} representing the on-shell particles. Note that the vertex function in Eqs. ([8]) 
and (19]) include the normalization factor of the wave function and momentum distribution 



function, given by |l9[ |: 



[zM)- (11) 



K -A;ion 72 Mo 

Note that in Eq. ( ITTj) . we have take = mg, i.e., m„ = for tt. To calculate the 
matrix element in Eq. ([9]), we choose a frame with the transverse momentum = so 
that = p^p^ > covers the entire range of the momentum transfers. Here, the relevant 
quark momentum variables are 

kt = {l-z')q+, kt = z'q^, k^^ = (1 - z')q^ + k'^, k2± = z' - k'^ . (12) 

By considering the good component as "/i = +", the hadronic matrix elements in Eq. ([2]) 
can be rewritten as: 

( 0|s7+75M|7r(p + g) ) = iU{p + qV , 
(7(g)|M7+75C?|7r(p + g)) = -e^^ (el-g±)p+, 

{^{q)\u^+d\nip + q)) = -^e^e^^e*g,p+ . (13) 

Using Eq. f lT2|) . the trace part in Eq. (fTOj) can be carried out. By comparing the last two 
equations in Eq. ( !T3l) with those in Eq. Q, we derive 

^ , 2^ . fdz'd^k^^, 1 flmd + E/c^e 2m„-AA;2e 



Fv{f) = -Am.j'-^^{z,kl)^_ 



— z 



J l md + {I- z){md-mu)klQ 2 m„ - z (m^ - m^) fc^G \ 
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where 



A = z{l- 2z'){rnd - rriu) - 2z'mu , 

B = z{l- 2z')md + md + (1 - 2z'){l - z)mu 



3 



e 



1 d^z, kl 



<l>{z,kl) dkl ' 



z = z' ( 1 - ^] , k = {k±, 4) , 

We note that to evaluate the form factors, we have to fix the meson scale parameter cu^r in 
the meson wave functions by fitting the meson decay constant, given by jsC 

rr- f dzd'^ki ^, , , m„ , 

W*'^- '=-'1(1^1) ■ 

B. Tensor Form Factor 

The tensor interaction is given by 

>Ct = ^ sin OJt ( u cr^,75 ) [ e a^'^ (1 - 75) ] • (17) 



The tensor form factor is defined by 14| : 



{l{q)\u(y^,y-fzd\'K{p^) ) = -«7rr (^^^'^ ~ ^t^^t) ■ (IS) 
For the LF good component of "/i = +" , one rewrites Eq. ( |T8|) as 

{l{q)\ua+'l,d\7r{p + q)) = (e*, ■ q^) . (19) 

At the quark level, the hadronic matrix element in Eq. (fT8|) is found to be 

{l{q)\u(rf,ul5d\7r{p + q)) = 



which leads to 



{l{q)\u(yp.l-,d\Tx{p\ q)) = 



7 



where 



/^Ifc-^ = Tr|75(-^2 + m„)2e,^*(^3 + m„)a^75(Vi + md)| . (22) 



From Eqs. fll9l) and fl^Tl) . we obtain 



h J 2 {2tt 



where 



Ci = -—. ^^7^ ^1 (1 - 2z' + 2z" - z'z){z' + z - 2z'z)kl 

z'z[i — zy[i — z') 1 

+ (1 - z)mu [2z'z{l - z')md + {z' + z + 2z''^z - Az z)md] 



(z - 



Z'z[l — -2) (1 — Z ) \ 

-{1- z){z' + z- z'z)miy (24) 



At the maximal recoil of = 0, we have 

^r(O) _ A f \2{l- z)k'^ + mu[zmd + {I- z)mu) 



h J 2{2txY ^ ' "^3 z{ml + k 

1 
3 



+ ^(m„ ^ rud) \ . (25) 



C. Numerical results 



To compute numerical values of the form factors in the LFQM, the uj parameter in the 
light-front wave function is fixed by other hadronic properties. For example, by using the 
decay constant of /tt = 130 MeV and the quark masses of rrtu = rrid = 250 MeV, we obtain 
cOtt = 301 MeV from Eq. (fT6|) . We note that this value of is just a typical one and its 
uncertainty mainly arises from those of the light quark masses. The transfer momentum 
dependences of Fy and Fa are shown in Figs. [T] and IH respectively. Note that the behaviors 
of the figure's sharps are independent of the quark masses but the values of Fy and Fa at 
the maximal recoil of = can be quite different as shown in Table [H The results in 
Table [T] also illustrate the main uncertainty from the quark masses. In Figs. [1] and [21 we 
have also included the experimental results fitted by the forms of Fv{p'^) = Fv{l + cxp^) 
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and Fa{p^) = Fa{0) with a constant parameter of a {4], while in Table. [U we list the values 
of FAy{0) in the ChPT and those from the data. We remark that the numerical values of 
the form factors at = for the pion case between the theoretical models seem to be less 
compatible comparing with those for the kaon in Ref. 10| in which the strange quark mass 
also enters. To illustrate the quark mass effects on the form factors in the LFQM, in Figs. [3] 
and m we plot three different sets of quark masses including the one in Figs. [1] and [2l It is 
clear that both Fv,a decrease as rriu^d increase. 
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FIG. 1. Fv{p^) as a function of the transfer momentum with niu = md = 250 MeV. 
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FIG. 2. Same as Fig. [2 but for Fy(p2 
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FIG. 3. Fvip^) in the LFQM, where ruu = ma = 230, 250 and 270 MeV correspond to sohd, 
long-dashed and short-dashed hnes, respectively. 




FIG. 4. Same as Fig. [3] but for Fv{p^). 

TABLE I. Values of ^4^(0) in (a) the ChPT at 0{p'^), (b) the ChPT at 0(/), (ci) the LFQM 
with i=l, 2 and 3 for rriu^d = 230, 250 and 270 MeV, repectively. 





(a) 


(b) 


(cl) 


(c2) 


(c3) 


Data [4] 


Fa(0) 


0.0112 


0.0102 


0.0151 


0.0131 


0.0113 


0.0117(17) 


Fv{0) 


0.0272 


0.0272 


0.02751 


0.0261 


0.0243 


0.0258(17) 
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The tensor form factor Eq. (123|) in the whole kinematic region for the LFQM is shown 
in Fig. 5. At = 0, we get the FT{0)/fT =0.220, 0.210 and 0.202 for m„,d = 230, 250 and 




FIG. 5. function of the transfer momentum p . 

270 MeV, respectively. 



III. DECAY BRANCHING RATIO 



In the TT"*" rest frame, we obtain the double differential decay rate as 



with 



^ |Mr = — fir(7r ^ ev)A, 



dxdy 2567r'^ ' ' 27r 
A = Aib{x, y) + Asd{x, y) + Aint{x, y) , 



Msix^y) 
AsD{.x,y) 



AiN{x,y) 



1- A 
<(1 - A) 



x^ + 2fl-rJ 1-x- 



-X 



Fy + Fa\ 



+ \Fv-FA\\y~X)] 



'Re[{Fy + FAr] (l-^-y) 



- Re[{Fv - Fa) 



A 



(26) 
(27) 



(28) 
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where x = 2£'^/m^, y = 2£'e/m,r, = ml/ ml and A = {x + y — 1 — re)/x. One can 
also relate the angle de^y between the e"*" and photon momenta with y and A. Explicitly, by 
neglecting the Ve, one has that 



A 



2 I de-y 



ysm \ — 



(29) 



The physical regions for x and y are given by 

<x< 1 



1-X + 



1 — X 



<y< 1 + Te 



(30) 



In Table II, we show the decay branching fractions of tt^ — )■ e^u^j in terms of the various 
contributions in Eq. for U = 130 MeV, m^ = 140 GeV I2I, m„,rf = 250 MeV and 

Br{n eu) = (1.23 ± 0.004) x 10"^ with the cuts of > 50 MeV and > 50 MeV. 
Note that in this kinematic region, the contribution from the SD part dominates the decay 
rate, which is sensitive to the V — A structure as well as new physics. We note that the 
total branching ratio in the LFQM are 2.937 and 2.320 x 10"^ for m^^d = 230 and 270 MeV, 
respectively. 

TABLE II. Decay branching ratio of vr — )• ez^e7 (in units of 10-S) in (a) the ChPT at O(p^), (b) the 
ChPT at 0(/) and (c) the LFQM with the cuts of > 50 MeV and > 50 MeV, respectively. 



Model 


IB 


SD 


INT 


Total 


(a) 


3.692 X 10^1 


2.356 


2.536 X 10-3 


2.727 


(b) 


3.692 X 10~i 


2.309 


2.850 X 10-3 


2.679 


(c) 


3.692 X 10"^ 


2.250 


1.840 X 10-3 


2.621 



In Table III, we give the decay branching ratio of vr — )■ ez/e7 in various kinematic energy 
regions in (a) the ChPT at 0(/), (b) the ChPT at 0(/), (c) the LFQM, (d) the green 
function method 16] and (e) the ChPT with a large Nc expansion [l7| as well as the data 
in Ref. j4|. Here, we have used ^ = 250 MeV in the LFQM. The errors in the parentheses 
of our results in Table III are from the decay of tt — eu. However, it should be noted that 
large uncertainties could arise from the various normalized coupling constants and the light 
quark masses in the ChPT and LFQM, respectively. In Fig. [6l we display the spectrum of 
the differential decay branching ratio as a function of x = 2ii^^/m^ in the ChPT at both 
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TABLE III. Decay branching ratio of tt — )• ev^l (in units of 10"®) in (a) the ChPT at 0{p'^), (b) the 
ChPT at 0(/), (ci) the LFQM with i=l, 2 and 3 for mu,d = 230, 250 and 270 MeV, repectively, 
(d) the green function method [3] and (e) the ChPT with a large Nc expansion as well as 
the data in Ref. [J| in various kinematic energy regions (in units of MeV) . 
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Data [4] 
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(cl) 
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(d) 


(e) 


50 


50 




2.614(21) 


2.727(9) 


2.679(9) 


2.85(8) 


2.62(8) 


2.29(8) 


2.81(38) 


2.58(8) 


10 


50 


40° 


14.46(22) 


15.04(5) 


14.99(5) 


14.93(37) 


14.63(37) 


14.19(37) 


15.08(58) 


14.77(40) 


50 


10 


40° 


37.69(46) 


38.28(13) 


38.12(12) 


38.29(12) 


37.87(12) 


37.37(12) 


38.4(10) 


38.9(9) 


0.5 


10 


40° 


73.86(54) 


76.66(25) 


76.31(25) 


73.67(22) 


73.57(22) 


72.58(22) 







0{p^) and 0{p^) and the LFQM. From Table III, we see that the results of the ChPT at 
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FIG. 6. Differential decay branching ratio as a function of x = 2E^/m.,^. 



0{p^) and 0{p^) are higher than those of the experimental data, which can be understood 
from Table I and Eq. (EHl) since the values of Fv-a{0) = Fv{0) - Fa{0) in the ChPT are 
larger than those fitted in the experimental data. Since the result from the LFQM agrees 
well with the data ^, it could lead to a strong constraint on new physics. We now examine 
the contribution to the decay from the tensor interaction in Eq. f[T71) . From Eqs. f|T7j) and 
(fTSl) . one obtains the new tensor contribution as 



Mr = te^ sin 6, {Ft e^q,) [ la^'{l - 75) ] 



(31) 
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Due to the above tensor interaction, Eq. (127|) should be rewritten as follows: 

A = AjBix, y) + Asd{x, y) + Aint{x, y) + Att{x, y) + Aibt{x, y) + Asdt{x, y) , (32) 
where the new terms are given by 



ATT{x,y) 



mz 



,2 



Ajbt{x, y) = 2^Re{FT) (l + r, - A - ^) 



A 



SDT 



m: 



^Re{FT){Fv - Fa) A(1 - A) . 



(33) 



Integrating over x and y variables in Eq.f lS^ and using the form factor in Eq.f E^ . we get 
the tensor related parts of the branching ratio as shown in Table. IV. To evaluate the 

TABLE IV. The tensor related parts of the decay branching ratio for tt — )• ev^r) (in units of 10^^) 
in the LFQM with the cuts of > 50 MeV and > 50 MeV. 



TT 


SDT 


IBT 


4.636 X 102 /2 


9.817 X 10-2 


2.536 X 10 /t 



total branching ratio including the tensor part, we can combine the results in Tables II and 
IV. By comparing the final result with the experimental data of — )■ e^z/e7, we extract 
/t = (3.48l|^^) X 10-^ shown in Table |V] in the LFQM for m^,^ = 250 MeV. In the table, we 

TABLE V. Form factor of fr in units of 10"^. 



LFQM 


[4] 


[m 


[13] 


[15] 


[16] 


3.48l|i 


-0.6 ±2.8 


-56± 17 


372 ± 120 


-115 ±33 


1 ± 14 



have also given other results in the literature including the single tensor form factor fitted 
by PIBETA [5|. We note that our resuh in the LFQM and that by PIBETA correspond to 
-1.0 X 10-3 < /t < 1.66 X 10-3 and -5.2 x 10"^ < /r < 4.0 x 10"^ at 90% C.L., respectively. 



IV. CONCLUSIONS 



We have studied the momentum dependent tt — 7 transition form factors F^yip^) in 
the ChPT and LFQM. In particular, we have found that ^^(0) = 0.0112, 0.0102, and 
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(0.151,0.0131,0.0113) and Fy(0) = 0.0272, 0.0272, and (0.0275,0.0261,0.0243) in (a) the 
ChPT 0(/), (b) the ChPT 0(/), and (c) the LFQM with m„,rf = (230,250,270) MeV, 
respectively, at the maximal recoil of = 0, Based on these form factors, we have calculated 
the decay branching ratio of tt — )■ eu^'j. Explicitly, we have obtained that in the SM with 
the cut of > nie and > 10 MeV with the relative angle O^^ > 40°, the decay branching 
ratio is 76.66 ± 0.25, 76.31 ± 0.25 and (73.67 ± 0.22, 73.57 ± 0.22, 72.58 ± 0.22) x 10"^ in 
(a), (b) and (c), respectively, while the experimental measurement is 73.86 x 10^^ by the 
PIBETA Collaboration. Since our results fit well with the data, we have also derived a 
constraint for the tensor interaction to be —1.0 x 10^"^ < /t < 1-66 x 10^'^ at 90% CL. in 
the LFQM . 
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